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Overlap fermions, which preserve exact chiral symmetry on the lattice, provide a powerful tool 
for investigating the topological structure of the vacuum. Applying this formulation to zero- 
temperature quenched SU(3) configurations generated by means of the Luscher-Weisz action, we 
define the topological charge density with and without UV filtering and study its properties by 
looking at the density profile and the two-point correlation function. We observe that the density 
possesses global sign coherent structures, which get increasingly tangled as more and more modes 
are included. This change of the structure is also detected by the increasing negative tail of the 
two-point function. We also study the inverse participation ratio of the eigenmodes and discuss 
their dimensionality. 
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1. Introduction 



Topological excitations, phenomenologically modeled as an interacting ensemble of instantons 
and anti-ins tantons, are believed to play a prominent role in the low-energy behavior of QCD, 
especially in realizing the axial U(l) anomaly (large r\' mass) and chiral symmetry breaking. It 
is thus of great interest to understand the topological structure of the QCD vacuum from first 
principles and to clarify whether instantons are indeed providing a realistic description. 

For this purpose, we employ the overlap fermion Dirac (Neuberger-Dirac) operator defined by 

D=-(l + ^£=), X=D W -?-, (1.1) 

where D\y is the Wilson-Dirac operator (r = 1.0). We set p = 1.4 as the optimal choice. Overlap 
fermions possess exact chiral symmetry on the lattice and provide n_ + n + exact zero modes, 
Dy^ = 0, with n_ (n + ) being the number of modes with negative (positive) chirality JsYn = 
~Yn (Y5Vn = +V / n h )- The index is given by Q = n_ — n + . The non-zero modes with eigenvalue 
X, Dyx = X\\fi, occur in complex conjugate pairs A and X* and satisfy (yIjYsVx) = 0- The 
topological charge density q(x), which satisfies the index theorem with Q = £ x g(x), is defined by 



(-v) = -Tr KM ^Dix.x)) 



(1.2) 



where 'Tr' is taken for color and spinor indices. The role of the low-lying Dirac eigenmodes for 
q(x) is exposed by applying the eigenmode expansion, where the cutoff A cut implies a kind of UV 
filtering, 

«u(*) = - I (1-7W*)' (L3) 

where ci (x) = (x)ys ( x ) is the local chirality of the mode with eigenvalue A. Note that the UV 
filtering maintains the index theorem independently of the cutoff such that Q = Ly^a^C*)' while 
controlling the UV fluctuation of the density. This is because the index is computed only from the 
zero modes which actually occur all with the same chirality. 



2. Observables and results 



In this presentation we study the topological charge density as well as its two-point correlation 
function. Furthermore we investigate the inverse participation ratio (IPR) of the eigenmodes. We 
use an ensemble of zero-temperature quenched configurations generated by means of the Luscher- 
Weisz gauge action (see Table [l] and ref. [[I]]). This action is suitable for topological studies since 
dislocations are greatly suppressed. 

Let us first look at the topological charge density. In Fig. |l], we show the density in a given 
time slice on a 12 3 24 lattice at j8 = 8.10. 1 The red (green) color surfaces are the isosurfaces 
corresponding to q(x) = ±0.0005. Periodic boundary conditions are imposed. In Fig. ^ we show 
the same plot for the zero-mode contribution alone and for the full density (no cutoff). It is apparent 

1 3D movies available on request. 
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Table 1: Simulation details: the Liischer-Weisz gauge action is used 



p 


a [fm] 


V=L 3 T (fm 4 ) 


# of conf . 


# of eigenmodes 


8.45 


0.095 


12 3 24 (3.38) 


116 


O(50) 


8.45 


0.095 


16 3 32 (10.6) 


267 


O(140) 


8.45 


0.095 


24 3 48 (54.0) 


186 


O(160) 


8.10 


0.125 


12 3 24 (10.1) 


254 


O(140) 




Figure 1: Distribution of the topological charge density in a given time slice of a 12 3 24 lattice at j3 = 8.10: 
cutoff Im A cut =0. 14, 0.28, 0.42, 0.56 (from left to right), isosurfaces at q(x) = ±0.0005 (red/green). 




Figure 2: The same plot as in Fig. 1 for the zero mode contribution (left) and the full density (right). 

that the topological charge density possesses global sign-coherent structures (clusters). The sizes 
of the clusters increase as the cutoff Ac Ut is relaxed. This feature is common also to the other 
configurations. We will address the local tomography of the clusters of topological charge in a 
future publication. 

Let us next discuss the two-point function of the topological charge density, 

C *(*) = 9{x)q{y)) , (2.1) 

' X 

where R = \x — y\ is the Euclidean distance (R = r/a). The result is shown in Fig. ||. The cutoffs 
chosen here are corresponding to those in Fig. [jj We observe that while there is always a positive 
core near the origin at any value of the cutoff, the tail is negative for sufficiently large cutoffs. As 
the cutoff increases, the peak of the positive core grows and the negative tail also increases. In 
fact, this is an expected behavior from the density distribution for various cutoffs as seen in Figs. |l] 
and ||. Since the density of both signs get increasingly tangled when more and more modes are 
included, any point of positive density comes closer to points of negative density. It is interesting 
to note that the topological susceptibility Xtop = (Q 2 ) /V is not affected by the cutoff. Therefore, as 
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Figure 3: Two-point function of the topological charge density qi on the 12 3 24 lattice at j8 = 8.10 defined 
for various A cu t, compared to that of the full density q (N con f = 3). The two-point function is averaged over 
pairs (x,y) with the same distance R. The right figure focusses at the region in R where the cutoff-related 
correlator turns negative. We obtained a similar picture at /3 = 8.45. 
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Figure 4: Two-point function of the full topological density on the 12 3 24 lattice at j3 =8.10 (iV C onf = 3) and 
on the 16 3 32 lattice at j3 = 8.45 (N con f = 2). The right figure shows the same as the left but magnified. Since 
only a few configurations were available for this comparison we did not estimate the errors. 



the core peak increases, the contribution from the negative region must grow to maintain the sum 
over C q (R) with respect to R. The dependence on the lattice spacing a is examined in Fig. |] for 
the full density, where the scale is introduced by multiplication with a -8 . Clearly, as a decreases, 
the core peak (size) is growing (shrinking), which seems to approach an expected behavior in the 
continuum theory although the finite positive core at our highest /3 stays at a radius R 2 a 
showing the limits of locality of the operator q(x) itself (see also ref. 

Finally let us discuss the inverse participation ratio (IPR) of the eigenmodes. The IPR of an 
eigenfunction (x) is defined by 

/ = ^I>aM, (2-2) 

xeV 

where Y,x€V Px( x ) = 1 w i tri Px{ x ) = Wx{ x )Wx{ x )- The IPR, / = \/f, denotes the inverse fraction 
/ of sites forming the support of Px(x) Limiting cases are / = 1 for p(x) = l/V (delo- 

calized), / = V for p(x) = 8 XrXQ (extremely localized). From the dependence on a and V one may 
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Figure 5: IPR: the dependence on ImX (left), the a dependence (right upper), and the V dependence (right 
lower). 

infer the effective dimension d of the eigenmode based on the relations I oc a d ~ 4 and I oc y 1 -^/ 4 . 
This means that there is no dependence on a and V if the mode is extending in d = 4 dimen- 
sions. In Fig. H (left) we show the IPR as a function of the (imaginary part of the) eigenvalue, 
where the IPR values are averaged in bins with bin size AImAj mp = 50 MeV, except for the zero 
modes which are considered separately. Here, Ai mp is the eigenvalue of the improved Dirac operator 
Amp = (1 — aD/2p)~ l D, which projects the eigenvalue of D stereographically onto the imaginary 
axis. The continuous eigenvalues of D[ mp still come in pairs ±/A, while the zero eigenvalues re- 
main zero. In Fig. ||, we also plot the averaged IPR with respect to its a dependence at fixed V 
(right upper), and with respect to the V dependence with a fixed a (right lower), separately for zero 
modes (7o) and lowest non-zero modes (7i ow ), the latter for < |Im Xi mp \ < 50 MeV. These plots 
exhibit that the zero and lowest-lying modes have an IPR clearly larger than one and are sensitive 
to the change of a and V, respectively. This indicates that these modes are localized to some extent 
and possess an effective dimension less than d = 4. From the V dependence, I = c\ +C2V l ~ d l A , we 
find that d = 2 and d = 3 describe the IPR data well for the zero and lowest-lying modes, respec- 
tively. On the other hand, for the higher modes (|Im Aj mp | > 200 MeV) we find I < 2, practically 
independent on a and V. Hence we may conclude that these modes are delocalized and extending 
in all d = 4 dimensions. The behavior of the IPR observed here is qualitatively the same as found 
in the SU(2) case [|]]. However, the actual values of the IPR are essentially smaller in our case and 
our fit results are also different. 



3. Summary 

Overlap fermions preserve exact chiral symmetry on the lattice and provide exact zero modes, 
which allow us to unambiguously compute the index Q of vacuum configurations. It is then possible 
to define the topological charge density q(x) and to expose the role of low-lying eigenmodes in q(x) 
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by applying the eigenmode expansion (UV filtering). This approach is gauge invariant and leaves 
the lattice scale unchanged in contrast to the cooling method. In this sense, overlap fermions are a 
powerful tool for investigating the topological structure of the vacuum. 

Using this formalism, we have investigated the localization properties of the topological charge 
density by looking at the density profile as well as the two-point correlation function on zero- 
temperature quenched configurations. We have also examined the IPR of the eigenmodes and their 
dimensionality. We have found that the topological charge density possesses global sign coherent 
structures, which get increasingly tangled as more and more modes are included. Such a change 
of density structures has also been detected in the behavior of the two-point function. We note that 
further preliminary studies show that the full topological charge density has a lower-dimensional 
laminar structure, together with a lumpy structure inside the sign coherent regions. Details will be 
presented in our forthcoming publication. The IPR has indicated that the zero and low-lying modes, 
typically for Im Aj mp < 200 MeV, are localized to some extent and have an effective dimension less 
than four. 

A similar investigation at finite temperature near T = T c on Nf = 2 dynamical gauge configu- 
rations is discussed in ref. 
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